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The asymptotic iteration method (AIM) is applied to obtain highly accurate eigenvalues of the 
radial Schrodinger equation with the singular potential V{r) = r^ + X/r" {a, A > 0) in arbitrary di- 
' mensions. Certain fundamental conditions for the application of AIM, such as a suitable asymptotic 

' form for the wave function, and the termination condition for the iteration process, are discussed. 

, Several suggestions are introduced to improve the rate of convergence and to stabilize the compu- 

■ tation. AIM offers a simple, accurate, and efficient method for the treatment of singular potentials 

' such as V{r) valid for all ranges of coupling A. 

o 

00 

O ■ PACS Nos.: 03.65.Ge. 



I. INTRODUCTION 

Attractive potentials with a strong repulsive core are common in atomic, nuclear and molecular physics pj — p^. 
A typical class of such potentials — IHl are the spiked harmonic oscillators V{r) = + The potential V{r) is 
called 'spiked' because of its pronounced peak near the origin. For a > 2, the potential is of relevance to quantum field 
theory, describing so called supersingular interactions for which matrix elements of the perturbation in the harmonic- 
oscillator basis diverge: thus every term in the perturbation series is infinite, and the perturbation expansion does 
not exist Q. There are several other reasons for interest in this class of potentials. First, it represents the simplest 
example of certain class of realistic interaction models in atomic, nuclear and molecular physics. Second, in the 
one-dimensional case, the perturbed-oscillator operator H — + V{r) — —cP/r"^ -f -I- A/r" — Hq + A/r", where 
p = —id/dr, may not converge to the original one Hq as A ^ (The Klauder phenomenon). Third, the perturbation 
series is ordered in fractional power [3, l^l in A. Since the early study of Harrell on singular perturbation theory for 
the ground-state energy of the Hamiltonian H, an enormous amount of work has been done to investigate the spectral 
problems of this operator. Most of the work, however, is either devoted to the study of the ground-state energy 
for A near zero, which represents the most challenging problem, or for particular values of the potential parameters. 
Owing to the difficulties inherited by the potential structure near the origin, much attention must be paid to the 
selected method for tackling these problems. The purpose of this letter is twofold. First to develop a simple and 



easily adopted technique, based on the asymptotic iteration method 2^1, to compute the eigenvalues of the radial 
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Schrodinger equation 



^^^^^ + 4 = Ei:{r),re [0,oo), 



(1) 



where a > 0, and the eigenfunctions satisfy the Dirichlet boundary condition ■0(0) = 0. This is valid for the 

one-dimensional case as well as for higher dimensions iV > 1 through 7 = Z + i {N — 3) , regardless of the values of 
the potential parameters. Second, to point out the importance of the correct form of asymptotic wave function to 
stabilize the iteration technique, and also to provide some suggestions to improve the rate of convergence of AIM 
when it is used to tackle Schrodinger eigenvalue problems with a wide variety of other singular potentials of physical 
and chemical interest. In order to achieve these goals, we develop first a wave function with the right exponential tail 
and which satisfies the Dirichlet boundary condition at the origin. This optimizes and stabalizes the use of AIM for 
computing the eigenvalues. The asymptotic iteration method (AIM) was original introduced [l^l to investigate the 
solutions of differential equations of the form 



y" = Ao(r)y' + so{r)y, 



(' = -) 



(2) 



where Ao(r) and so{r) are C°°-differentiable functions. Using AIM, the differential equation ^ has a general solution 



l{r) = exp - / p{t)dt 



C2 + C1 / exp / (Ao(t) + 2p{T))dT dt 



where, for sufficiently large n > 0, we obtain the p(r) values from 

Sn{r) _ Sn-l{r) ^ 

A„(r) A„_i(r) 



p{r) 



for 



An — -^'i-i + + AoA„_ 
Sn = s' 1 + soA„-i, n = 



1,2,3,. 



(3) 



(4) 



(5) 



It should be noted that one can start the iteration from n = with the initial condition A_i = 1 and s_i = 0. Since 
Ref. 22] the method has been adopted to investigate the spectrum of many problems in relativistic and non-relativistic 
quantum mechanics [l^] — [1^. In the process of applying AIM, especially in its application to eigenvalue problems of 
Schrodinger- type, such as ([1]), one usually faces the following two problems. 



A. Asymptotic solution Problem: 



The first problem we are confronted with in applying AIM is the conversion of the eigenvalue problem (the absence 
of first derivative) to standard form suitable to utilize AIM A general strategy to overcome this problem is 
to use an asymptotic solution ipaif) which satisfies the boundary conditions of the given eigenvalue equation. By 
substitution of the assumed exact solution with the form ip{r) — 4'a{f)f{i') into the eigenvalue problem, once the Aq 
and So have been determined, the sequences A„ and s„ can be computed using The energy eigenvalues are then 
obtained from the roots of the termination condition ([5]), which can be written in more convenient form as 

5n{r]E)^K{r;E)s,,^i{r;E)-K^i{r-E)s,,{r-E)=Q, n=l,2,.... (6) 

For Schrodinger's equation ([T]), the asymptotic solution is dominated by the harmonic oscillator wave function, since, 
for larger r, the dominant term of the potential is the harmonic oscillator term r^. The problem with such asymptotic 
solution, however, is that the behavior of the wave function near the origin has not been considered. 



B. Termination Condition Problem: 



The second problem results when the eigenvalue problem (now in the standard form for AIM application) fails to 
be exactly solvable. Indeed, if the eigenvalue problem has exact analytic solutions, the termination condition ^ leads 



Asymptotic Iteration Method for singular potentials 



3 



to an expression that depends only on the eigenvalues that is to say, independent of r. In such cases, the energy 
eigenvalues are simply the roots of 5n{E) = 0, n = 1, 2, . . .. For example, if a = 2, Eq.(IT|) is exactly solvable, and the 
termination condition ([6]) yields 



5n{E) =W{^i + ?, + 2-i - E), n = 0,l,2,3. 



(7) 



1=0 



The condition 5n{E) = leads to the exact solutions En = 4n + 3 + 27, n — 0, 1, 2, . . . as expected [23]- If the 
eigenvalue problem is not analytically solvable with the analytic form chosen, as for < a ^ 2, then the termination 
condition ([6]) produces for each iteration an expression that depends on both r and E. In such a case, one faces the 
problem of finding the best possible starting value r — that stabilizes the process so that it doesn't oscillate but 
converges. In principle, the computation of the roots of 5n{rvi\E) — should be independent of the choice of rg, 
nevertheless, the right choice of r, as we shall show in the present work, usually accelerates the rate of convergence to 
accurate eigenvalues E within a reasonable number of iterations. Generally, a suitable ro value is determined either 
as the location of the maximum of value of the asymptotic wave function, or as the position of the minimum value of 
the potential under consideration. A more general and systematic way to choose a suitable value for tq is still open 
question for further research. 

For the Schrodinger equation ([1]), we face these two problems. In the next section we develope an asymptotic wave 
function that satisfies the boundary conditions at zero and infinity. This asymptotic form is then used in section 3 
to initialized the asymptotic iteration method. In section 4, using a suitable value of ro, we exhibit and discuss the 
numerical results of AIM for a wide range of A and a > 0. Finally, in section 5, we comment on these results. 



II. ASYMPTOTIC WAVEFUNCTION FOR SINGULAR POTENTIALS 



For small r, one can neglect in ([T]) the energy E and the harmonic oscillator term as compared with the perturbative 
term A/r". Hence near the origin, ([1]) can be written as 



(Pi) 



(8) 



Using the transformation 

^(r) = V^0(<), t = /^r" 
a straightforward calculation shows that Eq.([8]) can be written as 



1 



a^(3 d<j) cr2^2 ^2^ 7(7+1) 



4J.3/2 

or, in more compact form, as 



-'^+3/2 dt 



d^ 
dt^ 



1 d(f> 



-2,7+3/2 



(27+1)' 



r3/2 



A 



„-Q + 2-2(T 



0. 



Therefore, for 2(T = 2 — a and [3 — VA/ct, wc obtain the modified Bessel's differential equation 



dt^ 



1 d(f> 
til 



f2 



1 



0, 



(9) 



where v = . The general solution of this differential equation is [3J] 

where ci and C2 are constants and Ii, and Ky are the modified Bessel functions of the first and second kind respectively 
[s^l- Consequently, an asymptotic solution of ([T]) for zero energy is given by 



^(r) = Vf[ci/,(/3r-) +C2X,(/3r-)]. 



(10) 
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We may now consider two cases, depending on the value of a. 

Case I (a > 2): In this case cr < 0, the boundary condition ijj{0) = forces ci = 0, hence 

2\/A 

M'r) = C2-V^K2:^{ -r^-*). (11) 

Q-2 a — 2 

From the asymptotic approximation 34] of K^{z), we know for large argument z that K^{z) = e^^/ \J^- Therefore, 
since 1 — ^ < 0, we have for small r, that 



Mr) = ci--J V °/^e =-2" . (12) 



Consequently, for a = 2m + 2, m > 0, 



i/;„(r) =r("+i)/2e-^, m>0 (13) 



up to a constant. In particular, if a = 4, (i.e. m = 1), one recovers the familiar limiting form of the solution when 
the repulsive potential is proportional to r~'^: 

y^a(r) = re >■ as r ^ 0. 
Case II (0 < a < 2): In this case cr > 0, the boundary condition i/'(0) — forces C2 = in (fTO| . thus 

2V^ 

ipair) = C2 ■ V?l2j+i{- r-i-^). (14) 

2-Q 2 — a 

For small argument z, the asymptotic approximation [34] of /^(z) = /r(!^ + 1) yields 

Mr) = r-'+^ (15) 

up to a constant. 

III. APPLICATIONS 

Depending on the degree of singularity of the potential at the origin, that is characterized by the positive parameter 
a, we have the following two cases: 

A. The case a = 2m + 2, m > 0: 

In this case the Hamiltonian operator 

H^-f,+r- + l^ + \, a>2 (16) 

leads to a non-Fuchsian singularity Q at r = of the Schrodinger equation Hip = Eil), because the potential term 
A/r" has a pole of order > 2. The asymptotic wave function developed in the previous section suggests that the exact 
solution of III), in the case a > 2, takes the form 

tpi^r) — e~^~'^^^ f{r). (17) 

The first exponential term in (jl7p takes into account the fact that for large r, the term in (jl6p dominates over all 
other terms of the potential, including, of course, 7(7 + l)/r^. For this wave function, Schrodinger's equation ^ now 
reads 

which is now amenable to AIM applications. Here, the primes of /(r) in (|18p denote the derivatives with respect to r. 
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B. The case a = 2m + 2, - 1 < m < 0: 
In this case, the exact solution of ([T|) assumes, for a < 2, the form 

i/.(r) =r'^+ie-'^/(r). 



(19) 



where again the exponential term in (flQl) takes into account that for large r, the term dominates over all other 
terms of the potential. In this case, Schrodinger's equation ^ reads 



fir) = 2 (^r - 2±1^ f\r) + [2^ + i + ^ - E ] f{r), 
which is suitable for an AIM application. 

IV. ITERATIVE SOLUTIONS 

For a given a, using (|18p or (pO|) . we can explicitly write Xo{r) and so{r) as: 
• a > 2: 

2\/A 1+ 



Ao(r) =(^2r-^- 
50 (r) = (^2 + m + ^ + (H2±i^I^^gl±l±!Ii) _ £; 



where m = ■^{a — 2). 
• a<2: 



Ao(r) = 2(^r-2±ij, 



(20) 



(21) 



(22) 



and, by means of the iteration formulas ([5]), we calculate A„(r) and Sn{i'), n = 1,2, . . .. The eigenvalues are then 
computed by means of the termination condition ([6]), namely 6n(r; E) — 0. With several symbolic mathematical 



TABLE I: The effect of using different values of r — on computing the eigenvalues using AIM for the Schrodinger equation 

~^d^ + ('"'^ + Tt) V' = An accurate value (to seven figures) is -Bcxact = 3.575 552. Here, N refers to the number of 

iterations needed to achieve such accuracy. 







N 


1 


2 


3 


4 


5 


15 


3.478854 


3.574644 


3.570009 


3.561477 


3.555018 


20 


Fails 


3.575335 


3.573786 


3.569455 


3.562441 


40 


Ull 


3.575551 


3.575505 


3.575298 


3.574726 


60 


Ult 


Fails 


3.575549 


3.575529 


3.575458 


85 


Ul") 


iLII 


3.575552 


3.575550 


3.575542 


90 


nil 


U11 


Done 


3.575551 


3.575545 


115 


Ult 


U11 


in; 


3.575552 


3.575551 


120 




iLII 


U5? 


Done 


3.575551 



programs available {Maple, Mathematica, etc) , the computation of the eigenvalues by means of the iteration method, 
provided it is set up correctly, is a straightforward calculation, even for the higher iteration steps. Most of our 
computations in the present work were done using Maple version 9 running on an IBM architecture personal computer 
(Dell Dimension 4400). As we mentioned above, the computation of the eigenvalues by means of (O should be 
independent of the choice of r. However, in some applications, for certain values of r, we may encounter oscillations 
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of the computed roots that seem to diverge in behavior. This is presumably due to rounding and computational 
errors in the algorithms used. In Table 1, we show the effects of choosing different starting values of r = tq on the 
number of iterations. It is clear that tq > 3 is sufficient as starting value of r. However, suitable choices for tq can 
significantly reduce the number of iterations needed to achieve the required accuracy. In many cases, we have removed 
the oscillating behavior by increasing the number of significant digits that Maple uses in numerical computations. In 
Table 2, we illustrate the effect of using different numbers of significant digits on the iteration convergence using a 
Maple environment. As indicated by the results in the table, a higher-precision environment can remove the oscillation 
behavior, as well as stabilize the numerical computation of the root problem by means of the iteration process. In 
order to accelerate the computation we have written our code for root-finding algorithm instead of using the default 
procedure Solve of Maple. For values of a = 1, AIM gives excellent results, even for extremely small value of the 

TABLE II: The effect of using different number of digits (in Maple 9) when computing the eigenvalues using AIM for the 
Hamiltonian — + + + Here we set r = 3 in ®. An accurate value is -Ecxact = 7.029816. Here A'^ refers to the 

ar 7 r ' — * 

number of iterations. 





Number of Digits 


N 


10 


14 


18 


22 


70 


7.0301659550 


7.0298160740 


7.0298162024 


7.0298162024 


75 


7.0298754200 


7.0298162383 


Done 


Done 


80 


7.0298152350 


7.0298162511 






85 


7.0298160500 


7.0298162636 




U11 


90 


7.0298162400 


Done 


a?; 


nil 



coupling parameter A. In table 3, we report the AIM results for a considerable range of A values. These computations 
have been made with — 3. For fractional a, such as 1/2,3/2, 1.9 etc., the method is still stabile and works well; 
however the number of iterations is much larger than that needed for integer a. The eigenvalues reported in Table 3 
are in excellent agreement with the exact eigenvalue computed by means of a numerical integration of Schrodinger's 
equation. 

TABLE III: The ground-state energy of Schrodinger's equation lU for a = 1 and different values of A using the present work. 
For these computations ro = 3. N is the number of iterations. 



a = 1 


A 


Eaim 


N 


1000 


190.723 307 439 784 825 395 54 


90 


100 


42.462 918 114 619 200 840 54 


32 


10 


10.577 483 539 371 157 357 99 


52 


1 


4.057 877 007 967 971 192 93 


64 


0.1 


3.112 066 906 502 466 751 74 


65 


0.01 


3.011 276 010 524 898 166 93 


62 


0.001 


3.001 128 301 284 079 220 13 


60 


0.0001 


3.000 112 837 137 807 781 38 


56 


0.00001 


3.000 Oil 283 783 881 865 84 


55 


0.000001 


3.000 001 128 379 089 204 58 


54 



In Table 4 we report the AIM results for the case a = 4, again for considerable range of A. The large number of 
iteration for small values of A reflect the stability of the method and also the applicability of AIM to treat such cases 
in one single formalism. Similar tables can be easily constructed for a = 3,5,6, etc. Fewer numbers of iterations are 
usually needed to achieve any required accuracy for the cases 7 > 0. For the ground state energy with 7 = 0, much 
attention to the value of ro must be paid to obtain accurate eigenvalues. An important observation is that: for large 
n, the computed roots by means of the termination condition Q are either in descending order or in ascending order; 
if an oscillation is observed, which appears to change the order, then the starting tq should be revised accordingly. 
The main point is that with the proper choice of rg, AIM is a stable and efficient method to obtain the eigenenergies 
to any degree of accuracy. 



Asymptotic Iteration Method for singular potentials 



7 



V. CONCLUSION 

The present work points out the importance of the asymptotic wave function used for initializing AIM sequences 
for Schrodinger eigenvalue problems. By introducing a wave function form that satisfies both boundary conditions at 
zero and at infinity, we were able to obtain accurate eigenvalues for the Schrodinger equation with singular potentials. 
Several suggestions are discussed to remove the numerical instabilities that may be encountered with direct utilization 
of AIM. Although we have focused our attention on the calculation of eigenenergies, the method also yields the 
corresponding eigenfunctions via equations Q and (|4|. 
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TABLE IV: The energy eigenvalues for Schrodinger's equation ([T]) with a = 4 and different values of 7 for a wide range of 
values of the coupling A. For these computations, ro varies over the range [4, 6.5] as A approaches zero. is the number of 
iterations. 



a = 4 


A 


7 




N 


1000 





21.369 462 532 163 464 497 98 


43 




1 


21.522 859 814 112 640 999 87 


39 




2 


21.827 883 093 646 909 321 84 


40 




3 


22.281 057 275 014 819 956 75 


38 




4 


22.877 334 674 778 463 023 90 


37 




5 


23.610 282 631 878 614 495 76 


36 


100 





11.265 080 431 752 838 088 14 


68 




1 


16.235 741 726 872 875 703 15 


69 




2 


16.801 763 365 978 787 670 47 


67 




3 


17.624 891 020 204 476 919 71 


65 




4 


18.677 251 190 728 170 156 30 


63 




5 


19.926 538 009 871 777 284 42 


58 


10 





6.606 622 512 024 943 661 69 


129 




1 


7.223 520 393 149 576 761 99 


124 




2 


8.352 483 528 249 905 501 41 


116 




3 


9.839 231 320 856 383 508 94 


100 




4 


11.544 000 451 519 493 138 54 


89 




5 


13.371 330 123 959 945 355 25 


78 


1 





4.494 177 983 369 188 


275 




1 


5.559 167 225 784 086 


246 




2 


7.224 287 163 959 573 


158 




3 


9.108 658 607 516 353 


131 




4 


11.062 241 719 384 166 


107 




5 


13.040 015 183 057 043 


88 


0.1 





3.575 551 992 


260 




1 


5.095 284 821 


200 




2 


7.025 961 149 


133 




3 


9.011 364 026 


90 




4 


11.006 336 099 


68 




5 


13.004 036 433 


50 


0.01 





3.205 067 495 


450 




1 


5.011 917 775 


356 




2 


7.002 658 316 


155 




3 


9.001 142 200 


85 




4 


11.000 634 789 


60 




5 


13.000 404 001 


51 


0.001 





3.068 765 


335 




1 


5.001 286 52 


259 




2 


7.000 266 58 


88 




3 


9.000 114 28 


53 




4 


11.000 063 49 


39 




5 


13.000 040 40 


33 



